Violation of critical universality at the antiferromagnetic phase transition of YbRh2Si2 
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We report on precise low-temperature specific-heat measurements, C(T), of YbRli2Si2 in the 
vicinity of the antiferromagnetic phase transition on a single crystal of superior quality (RRR~ 150). 
We observe a very sharp peak at Tjv = 72 mK with absolute values as high as C /T = 8 J/molK^. 
A detailed analysis of the critical exponent a around Tjv reveals a = 0.38 which differs significantly 
from those of the conventional universality classes in the Ginzburg-Landau theory, where a < 0.11. 
Thermal-expansion measurements corroborate this large positive critical exponent. These results 
provide insight into the nature of the critical magnetic fluctuations at a temperature-driven phase 
transition close to a quantum critical point. 



Thermodynamically driven phase transitions (PTs) oc- 
cur in manifold variants in solid-state physics. The Lan- 
dau theory of PTs was a major achievement in describ- 
ing PTs with one parameter, the order parameter [ij. 
Universal scaling dependences are valid for all continu- 
ous, i.e., second-order PTs and are only determined by 
the dimensionality of the system and the critical fluc- 
tuations In contrast, quantum PTs are driven by 
quantum fluctuations associated with Heisenberg's un- 
certainty principle 0, 0, 0| • At a quantum critical point 
(QCP), the ordering temperature of a continuous clas- 
sical PT is suppressed to zero temperature by a non- 
thermal control parameter (e.g., pressure or magnetic 
field), and the classical PT changes into a quantum one. 
Characteristics for a continuous quantum PT are strong 
deviations from Landau-Fermi-liquid theory in the phys- 
ical properties of metallic systems in a broad range of the 
phase diagram around the QCP. In heavy-fermion met- 
als, these non-Fermi-liquid phenomena have been stud- 
ied with great care, but a complete theoretical concept 
replacing the Fermi-liquidparadigm is still lacking (for 
recent reviews see Ref. 0, In this letter, we address 
the outstanding question whether a continuous magnetic 
PT driven by temperature in the vicinity of a QCP can be 
described within the framework of universality classes in 
the theory of PTs. To this end, detailed measurements of 
specific heat and thermal expansion of YbRh2Si2 around 
the antiferromagnetic (AFM) PT at Tn = 72 mK were 
performed. 

The heavy-fermion system YbRh2Si2 is well suited 
to study this interplay between quantum and classical 
PTs, because it is a clean, stoichiometric, and well- 
characterized metal situated on the magnetic side (T^v = 
72 mK), but extremely close to an AFM QCP, leading 
to pronounced non-Fermi-liquid behavior in transport 
and thermodynamic properties, such as the divergence 
of the electronic Sommerfeld coefficient 7 = C^^ /T, 
and a linear-in-T resistivity [3]. The observed temper- 
ature dependences disagree with the expectation for the 
three-dimensional (3D) spin-density-wave scenario that 



is applicable to many heavy-fermion materials at a QCP. 
Therefore, a new theoretical concept was developed, in- 
voking critical excitations that are inherently quantum 
[E 0|, explaining oj/T scalin g in CeCus.gAuo.i as well 
as diverging Griineisen ratios 13, 11 1 and a jump of the 
Fermi volume observed in Hall-effect measurements [l2j 
in YbRh2Si2. In this framework of 'local' quantum crit- 
icality, the Kondo effect is critically destroyed because 
local moments are coupled not only to the conduction 
electrons but also to the fluctuations of the other local 
moments Q . The recent observation of an additional en- 
ergy scale vanishing at the QCP which does neither cor- 
respond to the Neel temperature nor to the upper bound- 
ary of the Fermi-liquid region has boosted the interest in 
YbRh2Si2 [H. 

High-quality single crystals of YbRh2Si2 were grown 
from indium flux. The growing parameters were im- 
proved to gain larger single crystals with higher crys- 
tallinity. This new generation of single crystals exhibits 
residual resistivities as low as po ~ 0.5 ^Ocm, corre- 
sponding to a residual resistivity ratio of ~ 150 14]. The 
specific heat between 0.04 and 1 K was measured in a 
■^He/'*He dilution refrigerator, utilizing a semi-adiabatic 
heat-pulse method with background heating l^, 3| on 
a platelet-like single crystal (m ~ 10 mg) mounted on a 
silver platform. The addenda of the latter was measured 
separately and did not exceed 5% of the total measured 
specific heat at temperatures below 1 K. The determi- 
nation of the critical exponent necessitates a very sta- 
ble and accurate temperature of the environment. We 
achieved a temperature noise level below 30 /xK and a 
resolution of the temperature measurement better than 
±5/iK. For temperatures above 0.5 K, the specific heat 
was measured in a commercial PPMS with a •^He-insert. 
The 4/ contribution to the specific heat, C"^^ , was ob- 
tained by subtracting the non-magnetic, Clu, and the nu- 
clear, Cq, contribution from the measured specific heat, 
Cmeas, as C^^ = Cmeas - C Lu - Cq . C Lu was determined 
by measuring the specific heat of the non-magnetic ref- 
erence sample LuRh2Si2 below 10 K However, Clu 
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FIG. 1: (Color online) (a) 4/ increment to the specific heat 
plotted as C^f /T{T) on a logarithmic T scale for temper- 
atures above Tjv. Grey (red) line indicates a fit C/T oc 
ln(ro/r) with To = 30 K between 0.3 and 10 K. (b) C*^ /T{T) 
on a linear T-scale. The large peak at Tn = 72 mK shows the 
onset of AFM order. The peak sharpens dramatically upon 
reducing the power of the heat pulse during the measurement, 
leading from AT/T = 2 % (circles) to 0.1 % (squares), (c) Ab- 
solute values of C^^ /T at Tat as a function of AT/T. 



at 1 K contributes only to 1 % of Cmeas and the nuclear 
quadrupolar contribution, Cq = aq/T'^, calculated for 
Yb with aq — 5.68 x 10~^ JK/mol as determined from 
Mossbauer results [11 Ei, is well below 1% of C„,eas 
around Tn- 

In Fig. [T^, the 4/ increment to the specific heat is 
plotted as C^f /T on a logarithmic temperature scale 
above Tm- Below T - 10 K, C^^/T logarithmically in- 
creases with C^-^/T = a\n{To/T) with a = 0.15 J/mol 
and Tq — 30 K, in agreement with results on sam- 
ples of lower quality 0, Below T ~ 0.3 K, C^f /T 
increases stronger than logarithmically with a power law 
of C^/T cx T-°-34. This additional upturn was dis- 
cussed in detail in Ref. [l3| and ascribed to a 'break-up' 
of the heavy quasiparticles at the QCP. We note that this 
unique upturn is also visible in our new higher-quality 
sample, in contrast to measurements reported by Knebel 
et al. [2l|. 

Next, we focus on the specific heat at the magnetic PT. 
In Fig. [It, C^f/T is shown on a linear temperature scale. 
To determine the shape of the peak very accurately, we 
have reduced the power of the heat pulse to the low- 
est value at which an accurate determination of the heat 
capacity is still possible. In a standard experiment the 
power of the heat pulse is chosen such that the achieved 
temperature increase, AT, amounts to AT/T « 2%. 
The latter value was found out to be a good compro- 
mise of large signal-to-noise ratio and reasonable reso- 
lution The measurement with AT/T = 2% yields 
the black curve in Fig. [T]3. A pronounced peak is ob- 
served at Tjv. However, sharp features in the specific 
heat are artificially broadened when using such a rel- 



atively large AT/T. The measurements on the high- 
quality single crystal reveal a drastic sharpening of the 
peak when decreasing the power of the heat pulse such 
that AT/T — 0.1 %: The maximum value at T/v reaches 
/T « 8 J/mol K2. In Fig. [TJ:, this peak height 
is shown as function of AT/T. The absolute value of 
C"*//T at Tat increases strongly for AT/T < 2% but 
only slightly for AT/T < 0.5%, in contrast to what 
one would expect for a first-order PT. The scattering of 
the data below and above T^r is larger compared to the 
standard curve owing to the small temperature increase, 
AT < 0.1 mK. 

Well below T^, the Sommerfeld coefficient becomes 
constant, 70 = 1.7 J/mol K^, in agreement with previous 
studies |20| . However, the fact that C^f /T is larger below 
the magnetic PT than above is surprising and not yet 
understood. We note that the peak at T/v is remarkably 
sharp while, usually in heavy-fermion systems, anomalies 
in the specific heat associated with the onset of magnetic 
order are broad [H, |2§| . To the best of our knowledge 
there is no deep investigation of the critical exponent 
in specific-heat measurements on magnetically ordered 
heavy-fermion systems. 

In order to classify the magnetic PT, a detailed analysis 
of the critical behavior was done. We applied the usual 
fit function 
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to describe the critical behavior with the reduced tem- 
perature t = {T - Tn)/Tn; +{-) refers to t > {t < 0), 
respectively. The background contribution is approxi- 
mated by a linear t dependence {b + Et) close to T/v, 
while the power law (first term in Eq. [T]) represents the 
leading contribution to the singularity in C'^^{t) [13, H^- 
Often, a higher-order correction term is introduced in the 
fitting function; however, as the leading contribution in 
the case of YbRh2Si2 is very large, the data can be de- 
scribed well with Eq. ^ The procedure is as follows: 
The parameters a, , b, and E are fitted with fixed 
T/v = 72 mK, simultaneously for < > and t < 0. Then 
the fit is improved with an extended temperature range 
of typically 0.003 < \t\ < 0.1. This routine is repeated for 
slightly different T/v until the best fit (smallest deviation 
of the root-mean-square error) is obtained which, in ad- 
dition, gives a very accurate value of T/v- We performed 
this fitting procedure for three experimental data sets (1, 
2, 3) obtained with different AT/T (0.5%, 0.3%, 0.1%), 
the results of which are shown in Tab. HI For fit No. la 
we used the same experimental data set as for No. 1 but 
before fitting we subtracted the specific heat obtained 
at the critical field (T/v 0), in order to account for 
the quantum fluctuations. The critical field within the 
easy plane amounts to Be = 0.06 T and C^f /T grows as 
y-0.34 Y)q[o^ 0.3 K at Ba [2^. From Tab. HI it is obvious 
that the different experimental data sets lead to almost 
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TABLE I: Parameters derived from the fits of the specific- 
heat data for a high-quality single crystal of YbRh2Si2 around 
the AFM PT for different experimental data sets (see text). 
Fit fa is based on the same experimental data as fit f, but 
the contribution of the quantum fluctuations to the specific 
heat was subtracted prior to the fitting procedure. In fit 3, 
a Gaussian distribution of Tjv with STn /Tn ~ 3 ■ fO~* was 
included. 



TABLE IL Theoretically predicted critical exponent a of the 
specific heat for different universality classes (after Ref. ^2]). 
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FIG. 2: (Color online) Specific heat vs reduced temperature 
close to the AFM PT of YbRh2Si2. The data below and 
above Tn can be best fitted with a critical exponent of a = 
0.38 ±0.03. 



the same critical exponent, a = 0.38 ± 0.03. Moreover, 
a is not sensitive to the background subtraction which is 
completely different for fits 1 and la. 

The excellent quality of these fits becomes evident in 
Fig. [21 where the curve for fit 1 is shown below (solid 
line) and above (dashed line) Tn, together with the ex- 
perimental data (symbols). The total temperature win- 
dow in Fig. [2]corresponds to only 15 mK below and above 
Tn- {t) increases significantly close to Tm which is 
a first sign for a large positive exponent a. The fitting 
procedure described above yields a critical exponent of 
a = 0.38 ± 0.03, which describes the data in the entire 
temperature range around the PT. Also, no saturation 
(due to rounding effects of the PT) is observed at the 
lowest \t\ values. Only for 1ST jT = 0.1% we were able 
to measure below \t\ = where the data points tend 
to saturate. It may be attributed to residual sample in- 
homogeneities producing a microscopic smearing of Tm 
27|. We took these small rounding effects into account 
by a Gaussian distribution of Tm {STn/T^ = 3 • lO""*) 
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which has been carried out after the fitting procedure to 
Eq. ([1]) as described above. This method has been suc- 
cessfully applied to describe the rounding of specific-heat 
peaks [281 . [29|. The very narrow distribution ST^/Tn 
verifies the high quality of the YbRh2Si2 single crystal. 
The value of the critical parameter a is not sensitive to 
this additional parameter in fit 3, as can be seen from 
Tab.[Il 

In Tab. HI] we have summarized the theoretical values of 
a expected for the different universality classes, depend- 
ing on the symmetry of the system. It is obvious that 
the observed exponent strongly deviates from all possi- 
ble theoretical ones: The largest value for a is calculated 
for a 3D-Ising system, so that generally holds a < 0.11. 

To establish the anomalous critical exponent we have 
performed thermal-expansion measurements in the vicin- 
ity of Tn- The length change was registered perpendicu- 
lar to the tetragonal c axis on the same single crystal used 
for the specific-heat measurements, described above. For 
details of the measurement setup, cf. Ref. Thermo- 
dynamic relations reveal the same critical ex pon ent for 
the thermal expansion as for the specific heat [30| . How- 
ever, it is more accurate to extract the critical exponent 
by fitting a modified Eq. [T] directly to the measured rel- 
ative length change, AL/L, vs t Application of the 
same fitting procedure as for the specific heat gives a crit- 
ical exponent of the thermal expansion of a = 0.30±0.15, 
in good agreement with the specific- heat result. In Fig.[3| 
this fit is shown below (solid line) and above (dashed 
line) T/v, together with the experimental data (symbols). 
However, the minor variation of AL/L around Tjv (inset 
of Fig. [3]) severely complicates the analysis of the critical 
exponent and impedes a more accurate determination of 
a- The absolute temperature calibration is not as accu- 
rate as for the specific-heat measurements which leads to 
a slightly lower Tn — 68.5 mK. This, however, does not 
influence the analysis of the critical exponent, performed 
in the same way as in the case of the specific-heat data. 
Therefore, our thermal-expansion results do confirm the 
surprisingly large critical exponent in YbRh2Si2. 

The anomalous critical exponent a = 0.38 in YbRh2Si2 
raises the question if this critical behavior might be in- 
dicating a tendency towards weak first-order behavior. 
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FIG. 3: (Color online) Relative length change of YbRh2Si2 
as function of the reduced temperature close to the AFM 
PT. The data below and above Tjv can be best fitted with 
a critical exponent of a' — —0.7 ± 0.15 which results in a 
critical exponent of the thermal expansion of a = 0.3 ± 0.15 
[3H . The inset presents the same set of data on a linear t- 
scale. 

However, no first-order PT in YbRh2Si2 is reflected ei- 
ther in the imaginary part of susceptibility measurements 
on this new generation of single crystals, or by hysteresis 
in any transport and thermodynamic property. In addi- 
tion, no latent heat was resolved in our specific-heat mea- 
surements. Assuming a (somewhat smeared) first-order 
PT and integrating C^-f over a temperature window of 
\t\ < 0.1 one obtains H « ImJ/mol as an upper bound- 
ary of the latent heat which indeed would be difficult to 
resolve. 

Discarding of first-order PT, the value of a = 0.38 
leads to the conjecture that the AFM order of YbRh2Si2 
might belong to the new universality class of an XY an- 
tiferromagnet on a 3D stacked-triangular lattice, with 
a = 0.34 ± 0.06 [l^], as experimentally confirmed on 
CsMnBra [i^. However, for YbRh2Si2 such a scenario 
is more than unlikely, because this system has no ten- 
dency to chirality. The latter would require either a 
triangular arrangement of the magnetic ions or a non- 
centro-symmetric Yb site, allowing the Dzyaloshinsky- 
Moriya interaction to operate. As long as the magnetic- 
ordering vector of the magnetic phase and the disper- 
sion relation of the magnetic excitations are unknown, 
this question cannot be answered finally. Alternatively, 
the PT at Tn in YbRh2Si2 may well be influenced by 
the pronounced quantum fluctuations of the nearby lo- 
cal QCP which may substantially influence the spatial 
fluctuations of the classical order parameter and, thus, 
lead to internal degrees of freedom of the order parame- 
ter, preventing a description within the framework of the 
known universality classes. 

In conclusion, we have presented high-accuracy 
speciflc-heat measurements in the very close vicinity of 



the antiferromagnetic phase transition of YbRh2Si2 re- 
sulting in an extremely sharp peak at T/v = (72 ± 1) mK. 
A detailed analysis of the critical parameters yields an 
anomalous critical exponent a = 0.38±0.03 which cannot 
be explained within the conventional universality classes. 
Thermal-expansion measurements conflrm this violation 
of critical universality close to a continuous quantum 
phase transition. Therefore, the observed large critical 
exponent poses a considerable challenge for those theories 
describing the magnetic phase transition and its critical 
fluctuations in the vicinity of a QCP. 
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